Effect of superlattice modulation of electronic parameters on superconducting density 

of states in cuprate superconductors 
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Recent scanning tunneling microscopy on BSCCO 2212 has revealed a substantial spatial super- 
modulation of the energy gap in the superconducting state. We propose that this gap modulation 
is due to the superlattice modulations of the atoms in the structure, and hence the parameters in a 
microscopic model of the CuC>2 plane. The gap modulation is estimated using renormalized mean 
field theory for a t — t' — J model on a superlattice. The results compare well with experiment. 
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Scanning tunneling microscopy (STM) has become an 
important method toprobe the electronic structure in 
high-T c cuprates. [Hit S3 Recently STM has revealed 
a substantial modulation of the superconducting energy 
gap along the a-axis in BSCCO 2212 samples. [HQ The 
modulation appears in phase with the incommensurate 
periodic superlattice modulation which originates in the 
Bi-0 layers. @, 0] Structure studies show that this lat- 
tice modulation also affects the Cu02 planes, and so will 
modify the local electronic structure. Understanding this 
gap modulation may shed new light on our understanding 
of the superconductivity mechanism of high-T c cuprates. 
In this paper, we analyze the effect of modulating the 
parameters in a t — t' — J model on the superconducting 
state within the renormalized mean field theory (RMF). 
Q Our calculations show that the observed supercon- 
ducting gap modulation can be explained by a reasonable 
choice for the parameter modulation within a t — t' — J 
model as a consequence of the short coherence length, al- 
though the precise form of a priori t — t' — J modulation 
will need more concrete structure data. Zhu [l(| has pre- 
viously used a slave boson approximation to examine the 
effect of defect and acceptor induced parameter changes 
on superconductivity in the same model. Zhou et al. [11| 
used RMF keeping t, t' , and J fixed spatially but adding 
local Coulomb potentials to model dopant oxygen accep- 
tors close to the Cu02 planes to explain the nanoscale 
inhomogeneity found in STM experiments. In this paper 
we restrict our model to the superlattice structure. A 
somewhat different approach has been taken by Ander- 
sen, Hirschfeld, and Slezak [lj| who introduce a periodic 
modulation of the pairing interaction in a conventional 
d-w&ve BCS Hamiltonian to describe the superlattice. 

Shortly after the discovery of the BSCCO supercon- 
ductor, an incommensurate lattice modulation with a 
period of 5 unit cells along its a axis was found and 
characterized. @, [H Note this direction corresponds to a 
diagonal (1,1) direction in the Cu02 square lattice. The 
determination of the detailed displacements and site oc- 
cupation that occur in this superlattice modulation is dif- 
ficult to undertake in a complex material such as BSCCO. 



While there is unanimity among the various superlattice 
structures reported in the literature 
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substantial modulations occur in the key Cu02 planes, 
there is no agreement about their precise forms. [15| As 
a result, a priori estimates of the electronic parameter 
modulation within a t — t' — J model are still not possible 
at present. 

We start from a t — t' — J model on a square lattice in- 
cluding the no-double occupation constraint. Zhang and 
co-workers introduced the Gutzwiller approximation, re- 
placing the no-double occupation constraint by the clas- 
sical statistical weight factors g\ 3 and g\' 3 for the hop- 
ping and spin-spin superexchange processes, respectively. 
[9| The RMF that they derived reproduces the results of 
Variational Monte Carlo calculations on the full model 
quite well. Both have been proved to be capable of 
explaining many experiments qualitatively and in some 
cases even quantitatively [16l |. Decoupling the renormal- 
ized Hamiltonian gives a mean field form in the presence 
of pairing 
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where the operator c\ a (ci^ a ) creates (annihilates) an 
electron with spin a on the i th lattice site, tjj and t[ ■ are 
the hopping integrals for nearest neighboring (NN) sites 
(i,j) and next-nearest neighboring (NNN) sites 
respectively. Jy is the NN spin-spin coupling constant, 

/i is the chemical potential, and Xi,j = \ c \ a c i^ /> 



^i.j = ( c i.T c i l — c i.l c j-l) are t ne local NN particle-hole 
and particle-particle pairing fields. The Gutzwiller renor- 
malization factors gl' 3 and g l s ' J depend on the local doping 
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1 — m is the on-site hole concentration with 



where Si 
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As discussed above, the incommensurate lattice mod- 
ulation has a period of ~ 5 unit cells along the diagonal 
direction. In this paper we examine an electronic super- 
lattice structure with a commensurate periodicity, i.e., 
each unit cell consists of 10 sites as illustrated in Fig. 
[IJa). This superlattice has a period of 5y2 lattice con- 
stants along the (1,1) direction and y/2 perpendicular to 
it. The reciprocal lattice is shown in Fig. QJb). Wc 
use I (J) to index the supercell and i (j) for the origi- 
nal lattice. The lattice constant is set as the length unit 
with t = (0, ±1) or (±1, 0) for NN, and T = (±1, ±1) for 
NNN. Here we assume a simple cosine modulation, con- 
sistent with the dominant lattice modulation observed in 
Bi-O, Sr-, and Cu0 2 layers, 

hj+r/to = 1 + -[cos(Q-Ri + cos(Q-Ri +t )], (3) 

Ji,I+t/Jq = {tl,I+ T /to) i 

*i,i+r/*o = ! + tt [cos (Q • R I+r J + cos (Q • Ri+rJ] , 

where Q = (7r/5,7r/5), and Ri is the lattice coordi- 
nate within a unit cell. The average parameters are 
set to be t = 300 meV, Jo/to = 0.3, t' /t = -0.3, 
comparable to the experimentally observed values. A 
quadratic relation between t and J is preserved unless 
specified otherwise. The Hamiltonian in Eq. ([1]) with 
the supercell structure is solved self-consistently. The 
local density of states (LDOS) Ni (w) is determined by 
a Fourier transformation of the local Green's function 
iGi,a = (t t [c\ a (t) cl ct (0)]^ multiplied by a renormal- 

ization factor gl' 1 , where T T is the time ordering operator. 

Our simple model with modulations of t, t' , and J 
can reproduce many of the most prominent features of 
the STM experiments. The variation of the gap (about 
10%) within a wide doping range is reproduced by a rea- 
sonable choice of the electronic parameter modulation in 
our theory. It is worth noting that our model gives a 
very weak doping dependence of the relative gap vari- 
ance, which agrees well with the data recently reported 
by the Cornell group. [H, @] In addition, the low-energy 
subgap spectra are spatially homogeneous despite the in- 
homogeneity in the gap. Further the coherent peaks are 
always almost symmetrically located around zero bias. 
The gap is negatively correlated with the local doping 
concentration. In our calculation, we have normalized 
the gap value obtained from the RMF by a factor of 2 
to bring them in line with the more accurate Variational 
Monte Carlo results. Ildl 
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FIG. 1: (Color online) (a) The supercell structure. The black 
open dots present the sites in a unit cell with the coordinate 
inside the cell indicated, (b) The reciprocal lattice for a unit 
cell. 



To compare with the recent STM experiments, three 
typical doping levels with 5=0.13 (UD), 0.16 (OP), 0.19 
(OD) are considered. In our numerical calculation, a 
1000 x 400 supercell structure is introduced. The LDOS 
presented here is filtered by fast Fourier transformation 
with an energy window of 11 meV. We have studied the 
effects of various possible modulations of parameters t, 
t' and J to the LDOS and to the local gap values. The 
effects are qualitatively similar, although it is much less 
sensitive to vary t' than t or J. 

First we consider modulation of t and J while keeping 
t' unchanged, where the modulation amplitudes A = 4% 
and B = 0. The resulting variances of t and J are ~7 
and ~14%, respectively. The results shown in Fig. [2] are 
for doping 0.19, which are representative for all the three 
hole concentrations studied here. The LDOS is shown 
in Figs. HJa) and[^b). At low energy there is a clear 
homogeneous "V" shape [panel(c)] indicating the nodal 
structure (vf, va) is robust against this electronic mod- 
ulation. The homogeneity in LDOS at low energy was 
also found by Wang et ah, [13] in studies of disorder ef- 
fects. The coherent peaks with lower height are located 
at higher energy and the spectral weight suppressed at 
low energy is transferred to high energy near the band 
edge. Note that there is a multipeak character of the 
LDOS which is smeared by fast Fourier transformation. 
Similar multipeak character has also been observed in a 
recent STM experiment with high-energy resolution up 
to 2meV. [l8| The modulation of the superconducting en- 
ergy gap [shown in Fig. E^al)] has a variation about 10% 
for this parameter choice, which matches well with recent 
STM experiments. The doping variance achieved here is 
about 5%, i.e., 0.185-0.195 [shown in Fig. [^d)], compa- 
rable with the value of 7% recently reported in the STM 
experiments. 0, [lj)[ The negative correlation between 
the gap and the doping concentration is reproduced. In 
addition, there is a small deviation from the original d- 
wave symmetry due to the symmetry breaking imposed 
by the supercell structure orientation rotated by 45°. In 
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FIG. 2: (Color online) Panels (a)-(b) On-site LDOS with t 
and J modulated [A = 4%,B=0 in Eq. ©] with the peak 
positions indicated by arrows. The energy u is in unit of to 
(= 300 meV). Panel (c) The LDOS at low energies at two 
supercell sites (0,0) and (2,3), and panel (d) for on-site local 
hole concentration. The doping is 5 = 0.19. In (c)-(d), A — 
4% in the modulations of t and J modulation (red), and J 
only with t fixed (magenta), B = 45% in the modulation of 
t' (black), and A = 2% and B = 22.5% in the modulation of 
t, t' , and J (blue). 

Figs. G3al)-[2Icl) the modulations of the gap for different 
dopings are shown compared with the STM results. The 
phase <fi serving as the horizonal axis characterizes the 
electronic structure supermodulation which can be de- 
scribed by a consine form. The supermodulation phase 
mapping <j) (r) is consistent with the observed original to- 
pograph Z(r). @, @| The theoretical gap modulation vs 
Q • Rj (the range is linearly scaled to 360°) is shifted 
along the horizonal direction to match the gap minimum 
and maximum loci. Keeping the t and J modulation 
constant, a very weak doping dependence of the variance 
of the gap modulation is observed. Both the form and 
the magnitude of the gap modulation agree well with 
the experimental results on samples ranging between the 
underdoped and overdoped regime. The decreased agree- 
ment happened at low doping (5 = 0.13) is mainly due 
to the larger peak-to-peak distance of the on-site LDOS 
which possesses a multipeak structure before fast Fourier 
transformation. 

We have also considered the superlattice modulations 
of t' only, J only, and a combined modulation of t — t' — J. 
The low-energy LDOS and the local hole density for the 



FIG. 3: (Color online) The gap modulation due to electronic 
parameter modulations. The phase <f> characterizes the elec- 
tronic structure supermodulation which can be described by 
a cosine form. Experimentally the supermodulation phase 
mapping <b(r) is consistent with the original topograph Z(f). 
(Ref. [5j, |6() The theoretical gap modulation vs Q ■ R 7 (the 
range is linearly scaled to 360°) is shifted along the hori- 
zonal direction to match the gap minimum and maximum 
loci. The open diamonds with error bars in panels (al)-(cl) 
represent STM experimental data (Ref. [H, [(J) for three sam- 
ples (OD,OP,UD) with the energy gap A = 37, 47, 55 meV, 
respectively. Solid symbols are from our calculations, red cir- 
cles for t and J modulation (A — 4%, B = 0), black diamonds 
for t' modulation (A = 0, B = 45%), and blue squares for t, 
t' and J modulations (A = 2%, B = 22.5%). Color curves 
show analytic interpolations of our calculations between lat- 
tice sites. 

various choices of the modulated parameters are shown 
in Figs. Efc) andEJd), respectively. In Figs. I^a^c), 
we show the gap modulation obtained from the various 
choices of parameter modulations. It requires a substan- 
tial modulation of t', i.e., A = and B = 45%, corre- 
sponding to a range for t' = 0.165 — 0.435, to have simi- 
lar gap modulation with ~10% change in gap. The large 
variation in local hole density is due to the large modula- 
tion of t' required to obtain this modulation of the gap. 
By modulating J (A — 4%) alone with the quadratic rela- 
tion between t and J released, we have found only a slight 
difference from the case with both t and J modulated si- 
multaneously. Hence the results for modulating J alone 
are omitted in Fig. [3] due to the close similarity to mod- 
ulating t and J simultaneously. The low-energy behavior 
becomes more homogeneous and the variation of the dop- 
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ing concentration shrinks as shown in Figs. [2jc) and[2jd). 
Zhu who examined the nanoscale inhomogeneity [1 01 ] pre- 
sented a similar picture arguing that the presence of the 
randomly distributed out-of-plane dopant oxygens will 
lead to a change in the superexchange strength. From 
the cases studied above, we can reach a conclusion that, 
with similar variance of the superconducting energy gap, 
a modulation of t' leads to a larger modulation of dop- 
ing concentration while J and/or t tend to give a better 
low-energy homogeneity. We have examined a combined 
modulation for a choice of A = 2% and B — 22.5%. The 
LDOS, superconducting energy gap, doping concentra- 
tion, all show an additive behavior from the individual 
modulation of A and B with little interference, and the 
low-energy LDOS becomes more homogeneous. From our 
current study it is obvious that modulation of the elec- 
tronic parameters t, t', and J is able to generate the 
superlattice modulation phenomenon recently observed 
in STM. [|Q 

Pavarini et al. [2oT | have proposed that the varia- 
tion of the transition temperature T c between different 
cuprate families is controlled by the Cu-03 apical dis- 
tance through its effect on the NNN hopping parameter 
t'. The energy gap estimated by the RMF theory shows a 
similar dependence on t' , namely, the gap increases with 
increasing values of \t'\ at constant hole doping in agree- 
ment with the increase in T c argued by Pavarini et al. 
[20| It is therefore appealing to argue that this is the con- 
trolling feature of the superlattice modulation. However, 
there is a real problem with this interpretation. Slezak 
in his thesis || finds that if he uses the Cu-03 distance 



would be expected around a O 2 " ion but the enhanced 
gap suggests a locally reduced hole density. Several pro- 
posals have been made to reconcile theses two features 



I 111 Il9j . However, Yamamoto et al. [14j and also 
Eisaki et al. 22j conclude that in BSCCO a substantial 
number of Bi ions substitute for the Sr ions on planes ad- 
jacent to the CuOo planes. This suggests to us as it did 
to Eisaki et al. [22j, that since Bi 3+ acts as a donor ion, 
the regions of enhanced gap should be associated with 
Bi 3+ donors rather than O 2 " acceptors. 

We are very grateful to Seamus Davis and to S. Uchida 
for stimulating discussion. This work was supported by 
the RGC grant and Centre of Theoretical and Computa- 
tional Physics of HKS AR, RGC Central Allocation Grant 
of HKSAR, Visiting Professorship at The University of 
Hong Kong, and the MANEP program of the Swiss Na- 
tional Foundation. 



[1] 

[2] 
[3] 
[4] 

[5] 



modulation reported by Yamamoto et al. 14j there is a 
negative correlation, i.e. smaller values of the gap occur 
at larger Cu-03 distances. This is opposite to the conclu- 
sion of Pavarini et al. [2(| that \t'\ scales with the Cu-03 
distance. Thus the simplest interpretation that modu- 
lation of \t'\ is the controlling factor cannot be justified. 
Instead one cannot directly relate the reported structural 
data and the gap modulation. One can say, however, that 
the size of the modulations in the electronic parameters 
needed to explain the energy gap modulations are quite 
reasonable even if we cannot identify the specific values 
of individual parameter modulations. 

Lastly we comment on a couple of related issues. 
Slezak in his thesis [f| also reports on modulations of 
the bosonic mode fl accompanying the superlattice mod- 
ulations. However, as he points outs the isotope shift of 

has no effect on the gap which leads to the conclusion 
that both vary as a consequence of the structural modula- 
tion as was also suggested by Pilgram et al. [2l| A second 
feature of the BSCCO cuprates is the substantial disor- 
der in the STM gap maps. [3] This is characterized by 
regions with enhanced gap values associated with a high 
energy resonance in the STM spectra which has been in- 
terpreted as originating in a nearby dopant O 2 " ion. Q 
This leads to a problem since an enhanced hole density [16 



[10 
[11 

[12; 

[13: 

[14 

[15: 



S. H. Pan, J. P. O'Neal, R. L. Badzey, C. Chamon, H. 
Ding, J. R. Engelbrecht, Z. Wang, H. Eisaki, S. Uchida, 
A. K. Gupta, K.-W. Ng, E. W. Hudson, K. M. Lang, J. 

C. Davis, Nature 413, 282 (2001). 

K. M. Lang, V. Madhavan, J. E. Hoffman, E. W. Hudson, 
H. Eisaki, S. Uchida, J. C. Davis, Nature 415, 412 (2002). 
Michael Vershinin, Shashank Misra, S. Ono, Y. Abe, 
Yoichi Ando, and Ali Yazdani, Science 303, 1995 (2004). 
K. McElroy, Jinho Lee, J. A. Slezak, D.-H. Lee, H. Eisaki, 
S. Uchida, and J. C. Davis, Science 309, 1048 (2005). 
J. C. Davis, KITP conference: Strongly Cor- 
related Phases in Condensed Matter and De- 
generate Atomic Systems (May, 2007), website 
http://online.kitp.ucsb.edu/online/coldatoms07/davisl/ 
J. Slezak, Ph.D. Thesis, Cornell University, available at 
http: //people. ccmr.cornell.edu/~jcdavis / theses /james / . 
J. A. Slezak et al., (to be published). 
Yan Gao, Peter Lee, Philip coppens, M. A. Subramanian, 
A. W. Sleight, Science 241, 954 (1988). 

E. A. Hewat, M. Dupuy, P. Bordet, J. J. Capponi, C. 
Chaillout, J. L. Hodeau, and M. Marezio, Nature (Lon- 
don) 333, 53 (1988). 

F. C. Zhang, C. Gros, T. M. Rice, and H. Shiba, Super- 
con. Sci. Tech nol. 1, 36 (1988). 

Jian-Xin Zhu, cond- mat/0508646 (unpublished). 

Sen Zhou, Hong Ding, and Ziqiang Wang, Phys. Rev. 

Lett. 98, 076401 (2007). 

Brian M. Andersen, P. J. Hirschfeld, James A. Slezak, 
Phys. Rev. B 76, 020507 (2007). 

D. Grebille, H. Leligny, A. Ruyter, Ph. Labbe and B. 
Raveau, Acta Cryst. B 52, 628 (1996). 

Akiji Yamamoto, Mitsuko Onoda, Eiji Takayama- 
Muromachi, Fujio Izumi, Toru Ishigaki and Hajime 
Asano, Phys. Rev. B 42, 4228 (1990). 
For example, Grebille at al. [131 report a large vacancy 
wave on the 03 apical oxygen site. Yamamoto et al. jl4l ] 
find full occupancy on this site but substantial numbers 
of Bi and Ca-ions on the Sr-sites which are in the same 
plane. 

P. W. Anderson, P. A. Lee, M. Randeria, T. M. Rice, N. 



■5 



Trivedi, and F. C. Zhan, J Phys. Condens. Matter 16, 
R755-R769 (2004); Kai-Yu Yang, C. T. Shih, C. P. Chou, 
S. M. Huang, T. K. Lee, T. Xiang, and F. C. Zhang, Phys. 
Rev. B 73, 224513 (2006); Arun Paramckanti, Mohit 
Randeria, and Nandini Trivedi, Phys. Rev. B 70, 054504 
(2004). 

[17] Qiang-Hua Wang, Z. D. Wang, Yan Chen, and F. C. 

Zhang,, Phys. Rev. B 73, 092507 (2006). 
[18] A. C. Fang, L. Capriotti, D. J. Scalapino, S. A. Kivelson, 

N. Kaneko, M. Greven, and A. Kapitulnik, Phys. Rev. 

Lett. 96, 017007 (2006). 



[19] Tamara S. Nunner, Brian M. Andersen, Ashot Melikyan, 
and P. J. Hirschfeld, Phys. Rev. Lett. 95, 177003 (2005). 

[20] E. Pavarini, I. Dasgupta, T. Saha-Dasgupta, O. Jepsen, 
and O. K. Andersen, Phys. Rev. Lett. 87, 047003 (2001). 

[21] S. Pilgram, T. M. Rice, M. Sigrist, Phys. Rev. Lett. 97, 
117003 (2006). 

[22] H. Eisaki, N. Kaneko, D. L. Feng, A. Damascelli, P. K. 
Mang, K. M. Shen, Z.-X. Shen, and M. Greven, Phys. 
Rev. B 69, 064512 (2004). 



